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A STUDY OF BERYLLIUM AND BERYLLIUM -LITHIUM 
COMPLEXES IN SINGLE -CRYSTAL SILICON* 

By Roger K. Crouch, James B. Robertson, 
and T. E. Gilmer, Jr.** 

Langley Research Center 

SUMMARY 

When beryllium is thermally diffused into silicon, it gives rise to acceptor levels 
191 meV and 145 meV above the valence band. Quenching and annealing studies indicate 
that the 145-meV level is due to a more complex beryllium configuration than the 
191-meV level. When lithium is thermally diffused into a beryllium -doped silicon sample, 
it produces two new acceptor levels at 106 meV and 81 meV. Quenching and annealing 
studies indicate that these new levels are due to lithium forming a complex with the 
defects responsible for the 191-meV and 145-meV beryllium levels, respectively. Elec- 
trical measurements imply that the lithium impurity ions are physically close to the 
beryllium impurity atoms. The ground state of the 106-meV beryllium -lithium level is 
split into two levels, presumably by internal strains. Tentative models are proposed to 
explain these results. 


INTRODUCTION 

Technological advancements in recent years have led to an ever increasing demand 
for new and better solid-state electronic devices. To assist in the development of such 
devices, numerous studies on the electrical (refs. 1 to 4) and optical (refs. 5 to 8) proper- 
ties of the basic semiconducting materials have been carried out. As the technology 
progressed in growing and purifying single -crystal silicon, it became possible to intro- 
duce various concentrations of specific impurities into the crystal and then to catalog 
the properties of the silicon containing the known impurity. Initially, it was found that 
Group V impurities produced n-type silicon and Group HI impurities produced p-type 
silicon. It was shown that effective -mass theory (refs. 9 to 12) gave excellent agreement 
between predicted and measured values for the excited-state spectra of shallow donors, 

* The information presented herein is based on a thesis submitted by Roger K. 
Crouch in partial fulfillment of the requirements for the degree of Doctor of Philosophy 
in Physics, Virginia Polytechnic Institute and State University, Blacksburg, Virginia, 

1971. 

**Professor of Physics, VPI and S.U., Blacksburg, Va. 




that is, Group V impurities. (Shallow impurity states, either donor or acceptor, are 
those bound states whose binding energy is much less than band-gap energy.) The rather 
large discrepancies between predicted and observed ground-state energies were under- 
standable due to the simplified potential assumed in effective -mass theory (refs. 12 
and 13). When effective -mass theory was applied to shallow acceptors, it gave relatively 
good qualitative predictions for the excited states but it did not give good quantitative 
values. Acceptor states are quite complicated, and the theory is still undergoing devel- 
opment (refs. 14 and 15). 

Other studies showed that oxygen was present in large quantities in Czochralski- 
grown silicon; but the oxygen was electrically inactive in most cases and only modified 
the infrared -absorption spectra of silicon by introducing vibrational absorptions, pri- 
marily around 9 jum (refs. 16 and 17). The oxygen almost certainly occupied an inter- 
stitial lattice site in one of the six equivalent sites around a silicon-silicon bond (ref. 18). 

When lithium is introduced as an impurity in silicon, it forms a complex with the 
oxygen which alters the electrical and optical properties (refs. 19 to 21). The lithium 
occupies a tetrahedral, interstitial position and is described by effective -mass theory 
(ref. 21). In float -zone silicon, oxygen content less than 10^® cm - ^ the lithium gives 
rise to an energy level 32 meV below the conduction band. In Czochralski silicon, oxygen 
content typically 10*® cm"^, the lithium energy level is located 39 meV below the conduc- 
tion band. Further studies indicate that the 9-pm vibrational absorption of the inter- 
stitial oxygen is shifted to about 10 jum (ref. 22). Lithium has also been used to compen- 
sate heavily doped p-type silicon in order to reduce free -carrier absorption and to study 
pair formation between lithium and boron. Vibrational absorptions have been found which 
have an isotopic dependency on the lithium and the boron, and a tentative model has been 
presented for the lithium -boron complex (refs. 23 to 25). Some recent work deals with 
vibrational absorptions due to carbon and carbon-oxygen impurities in silicon (ref. 26). 

Group II impurities such as zinc and beryllium have been studied in germanium 
(refs. 27 to 29) and have been found to produce acceptor levels which, to a first approxi- 
mation, can be described by a helium model. (A helium model is essentially a system 
which can be treated as a doubly charged nucleus with two carriers (either holes or 
electrons depending on the nuclear charge) bound to it at low temperatures. The first 
electron (hole) will have excited-state energies comparable to those predicted by the 
hydrogen model; the se cond ele ctron (hole) will have a ground-state energy about four 
times deeper than the first electron. Also, the energy spacings between the excited 
states of the second electron will be larger by a factor of four.) The Group II impurities, 
magnesium and beryllium, have been studied in silicon (refs. 30 to 34). It is interesting 
to note that the magnesium produces n-type silicon and occupies an interstitial lattice 
site, while beryllium produces p-type silicon and is therefore believed to occupy a 
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substitutional lattice site. Magnesium introduces two related donor states described by 
a helium model. Beryllium introduces two independent series of infrared absorptions, 
which are separated by about 40 meV, neither of which seem to have the associated, 
deeper level expected for the helium model. (While it is uncommon for an impurity to 
produce two independent levels, beryllium is not unique in this respect since sulfur intro- 
duces at least two independent donor levels. See ref. 35.) When lithium is introduced 
into a beryllium -doped silicon sample, two new series of infrared absorptions are found 
at energies closer to the valence -band edge, indicating lithium forms a complex with the 
beryllium. 

A study of the electrical and optical properties of beryllium and beryllium -lithium 
complexes in silicon was done to better understand the unusual properties associated 
with these impurities. The results of this study are presented and discussed in the pres- 
ent paper. Tentative models for the two beryllium and the two beryllium -lithium impur- 
ity sites are proposed. 
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SYMBOLS 


numerical constants (eq. (6)) 

a 2 a 2 a 2 

Laplacian operator, — - + — - + — ^ in rectangular Cartesian coordinates 
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electronic charge 


energy level of shallow donor 


Ea 

Fj(F) 


acceptor ionization energy 
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hydrogenlike envelope wave function 

Planck constant 
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light intensity 


incident light intensity 


total angular momentum quantum number 


wave vector 


3 



m 


m 0 


n 

r 

R 


Rh 

t 

T 

u(D 

v(r) 


a 


“i 


x 



p 

T 


effective mass of charge carrier 
rest mass of electron 
longitudinal effective mass 
transverse effective mass 
free -carrier concentration 
radial vector 

reflectivity of sample surface ; 

Hall coefficient 
sample thickness 
temperature 

Coulomb potential due to interaction between electron and impurity ion 
periodic potential of an electron in a perfect silicon lattice 
absorption coefficient (eq. (7)) 

weighting factor for relative contribution of each equivalent minimum 
parameter (eq. (6)) 
dielectric constant 

spin-orbit splitting of valence bands at k = 0 

Hall mobility 

resistivity 

percent transmission 
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Bloch wave function at the jth conduction band minimum 


Vfr) 

i//(F) donor -state wave function 

i/a(F) impurity wave function 

SHALLOW IMPURITY STATES IN SILICON 
Donors 

When Group V elements are placed into a substitutional silicon lattice site, four of 
the electrons in the outer shell of the impurity atom are needed to complete the covalent 
bonds with the four adjacent silicon atoms. The fifth electron, however, is bound by a 
Coulomb attraction to the remaining charge on the impurity ion. If this system is treated 
as a hydrogen atom, where the potential energy of the electron is modified by the dielec- 
tric constant of the material and effective mass replaces the free -electron mass, then 
the excited-state spectrum associated with the impurity should be observed in the far 
infrared. Experimental data on the Group V impurities, phosphorus, antimony, and 
arsenic, convincingly demonstrate the existence of bound states as demanded by the 
hydrogen model for donors. 

The Schrodinger equation for the donor-state wave functions has the form 

^(r) = E <//(?) (1) 

where v(r) is the periodic potential of an electron in a perfect silicon lattice and u(F) 
is the additional potential due to the impurity ion. For large r 

U(r) = - (2) 

nr 

where k ~ 12 is the dielectric constant of silicon. In silicon, the energy -band structure 
of the conduction band has been found to have six equivalent minimums located along the 
<100> directions in k space (ref. 12). It has been shown by Kohn and Luttinger (ref. 9) 
that the impurity wave functions have the approximate form 

6 

^(r) = ^ o'.Fj(r)<p.(r) (i = 1, . . ., 6) (3) 

j=l 

where <^(F) is the Bloch wave function at the jth conduction band minimum and or j 
are the numerical coefficients which describe the relative contributions of each of the six 
equivalent minimums. The Fj(F) is a hydrogenlike envelope wave function which, for a 


_ *1 v 2 + y(r ) + U(r) 
2m 
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band minimum in the <001> direction, satisfies the following effective -mass wave 
equation: 



2m t 



d * 

2m 7 9z 2 
] 


e 

KT 


Fj(r)= EFj(r) 


(4) 


where m^. = 0.19mQ and m^ = 0.98mQ (ref. 36). 

When various Group V impurities are observed experimentally, the excited-state 
structure observed in the spectra are in excellent agreement with theoretical predictions 
resulting from approximate solutions of equations (4). Since the orbits of these excited 
states are quite large, the potential in equation (2) should describe the situation very well. 
However, the value of the ground-state energy is different for different impurities; but 
this is not unexpected, since for small values of r the assumption of a screened ion 
potential is very poor. Local strains, a varying dielectric constant, and numerous other 
complications arise in the immediate vicinity of the impurity ion. There are no adequate 
theories to predict the ground -state energies of the various impurities. 


Acceptors 


When Group HI impurities are placed on substitutional lattice sites in silicon, the 
unsatisfied covalent bond associated with one of the adjacent silicon atoms is an available 
electron state which can accept an electron from the valence band. When an electron is 
taken from the valence band, the resulting positively charged hole will have a Coulomb 
attraction for the negatively charged impurity ion. Thus, an effective -mass hydrogen 
model can also be used to describe the acceptor states. While these states are analogous 
to the donor states in a general way, the theory must incorporate the complexity of the 
valence band near its maximum at k = 0 (see fig. 1). 


If spin is not included, the band is threefold degenerate at k = 0; but when spin is 
included, the resulting sixfold degeneracy is split by the spin-orbit interaction into a 
fourfold degenerate P3/2 band and a twofold degenerate Pjy2 band. The P3/2 band 
corresponds to the atomic J = 3/2 states, and the P-^ band corresponds to atomic 
J = 1/2 states. 


Because of the nature of this degeneracy (i.e., the fact that three different bands 
are involved), the effective -mass theory is more complicated for acceptors than for 
donors in silicon. Kohn (ref. 12) has treated the problem for silicon and found that the 
total wave functions of the acceptor states have structures given approximately by 


6 

^ k (r ) = ^ Fj^rOcp-Cr) (5) 

j=l 
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Reduced wave vector 


Figure 1.- The spin-orbit split valence band of silicon with associated acceptor states; 
X is the spin-orbit splitting at k = 0; X = kb meV. 

where F^(r ) are hydrogenlike envelope functions modifying the Bloch functions cp .{ r ) 
The index k distinguishes the degenerate functions among the impurity states, and the 
functions Fj(r ) are determined by the following six coupled hydrogenlike effective - 
mass wave functions: 



+ fe F / r M E+Xe i) F i (F) 


3 = 1 , 2 ,. . .,6 

«,£ = 1,2,3 , 


( 6 ) 


where D..f are numerical constants characteristic of the material and have dimensions 

of -ft 2 /m; 6^ = 0 for j = 1 to 4 and = 1 for ] = 5 to 6. It may be shown that if X 

were very large compared to E, the six equations would separate into four equations 
associated solely with the p 3 / 2 states and two equations associated solely with the 
Pjy 2 states. In germanium, the spin-orbit splitting is about 30 times larger than the 
ionization energy, and only four equations are required to describe the acceptor states. 
For silicon, the energy of the acceptor states is on the order of 40 to 200 meV and X is 
on the order of 40 meV; therefore, the wave functions of all six valence bands will enter 
with appreciable amplitude into the total wave function of the acceptor states in silicon. 
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When the experimental infrared-absorption data are compared to the theoretical 
predictions, the agreement is qualitatively good; but there are small differences in the 
spacings of the excited-state spectra for different impurities (fig. 2(a)). Also, as in the 
case of donors, the assumption of constant k is invalid for small values of r; there- 
fore, the depth of the ground state is chemically dependent. The acceptor -state equations 
for silicon, as is also the case for germanium, have been solved only approximately by 
variational techniques. However, it turns out that the theory and experiment do not agree 
as well for acceptor states in silicon as for these states in germanium, presumably due 
to the complexity introduced by the small value of the spin-orbit splitting in silicon. 

As previously mentioned, the top of the Pg/2 va l ence band is fourfold degenerate. 
Many of the bound states associated with it are also fourfold degenerate. When external 
stresses are applied to silicon, the degeneracy of the localized states may be partially 
lifted. If the stress is applied uniaxially along the major axes of the crystals, it is theo- 
retically possible to determine the degeneracies of the ground states and the excited 



30 33 4 0 65 70 

Energy, meV 


75 


80 


(a) Lines associated with (b) Lines associated with 

shallow acceptors (ref. 8). beryllium complexes. 

Figure 2.- Absorption spectra of Group III impurities and the beryllium associated 
impurities. Vertical lines indicate the predicted location of lines using 
effective-mass theory; all ordinate values are schematic representations of 
relative absorption intensities. 
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states. Onton, Fisher, and Ramdas (ref. 8) have studied the Group III impurities (boron, 
aluminum, gallium, and indium) in silicon by using these techniques. They have shown 
that all of these impurities are located in a tetrahedral lattice position. The ground state 
associated with this position is fourfold degenerate; and when stresses of sufficient mag- 
nitude are applied, the level is split into two twofold degenerate levels. 

This section has presented a brief review of the theory for shallow impurity states 
in silicon. The ionization energies associated with the impurities reported here are large 
enough compared to the band gap energy to cause doubt as to the applicability of shallow - 
acceptor theory; however, the excited -state spectra of all the impurities studied appear 
to be in excellent agreement with the Group III impurities (see fig. 2). Because of this, 
at least the beryllium -associated impurity excited states have been regarded as shallow 
impurity states. 


DESCRIPTION OF THE EXPERIMENT 
Optical Studies 

Both Czochralski and float -zone single -crystal silicon (ref. 36) were used in pre- 
paring the test samples. Table 1 lists each sample reported, along with its method of 
growth, predominant impurity, resistivity prior to and after doping, doping time and tem- 
perature. The samples were polished disks 25.4 mm in diameter, and the thickness 
varied from about 0.5 mm up to 3.0 mm. Infrared -absorption spectra were obtained 
over a range from 1.0 to 50 pm using two grating spectrophotometers. Both instruments 
were purged by dry air during all runs (ref. 37). 

TABLE 1.- GENERAL INFORMATION ABOUT SAMPLES 


Sample 

identification 

Data 

presented 
in figure 

Method of 
growth 

(a) 

Major 
impurity 
prior to 
doping 

(a) 

Resistivity 
prior to 
beryllium 
doping, 
fl-cm 
(a) 

Beryllium doping 

Resistivity 

after 

beryllium 

doping, 

f2-cm 

(a) 

Time, 

hr 

Temperature, 

°C 

269-4 

5,6 

Float -zone 

Boron 

100 

i 

1300 

0.3 

X-9 

7, 8(b), 9 

Float -zone 

Phosphorus 

100 

i 

1300 

.25 

858-13 

8(a) 

Float -zone 

Phosphorus 

800 

i 

1300 

— 

Ox-11 

10 

Czochralski 

Boron 

140 

2 

1320 

— 

X-10 

11, 12 

Float-zone 

Phosphorus 

100 

1 

1300 

.2 

Ox-4 

13 

Czochralski 

Boron 

140 

2 

1340 

.4 

422-14 

14 

Float -zone 

Boron 

260 

ll 

A 

1340 

.2 

957-23 

15 

Float -zone 

Boron 

380 

1 

1300 

.5 


a Information supplied by manufacturer. 
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The infrared spectroscopy was performed at temperatures less than 20 K for all 
figures presented in this report unless otherwise noted in the figure. The low tempera- 
ture of the samples was attained through the use of a liquid helium Dewar in which the 
sample was clamped into a copper sample holder that was in contact with the cryogenic 
liquid. The temperature of the samples was not controlled in most cases and varied from 
around 8 K to 20 K depending upon the impurity concentration, thickness of the sample, 
and grade of thermal contact. Sample temperatures were measured with a commercially 
available germanium resistance thermometer. A piece of germanium was sealed inside a 
small copper block with a helium transfer gas inside. Four wires were mounted on the 
germanium to measure the resistance which had been previously calibrated as a function 
of temperature. The copper block was mounted against the silicon by using vacuum 
grease. 

All plots of optical data are presented in terms of percent transmission, rather than 
absorption coefficient, as a function of wave number. Since the absorptions are very 
strong in most cases, errors in percent transmission (up to 5 units) would lead to large 
errors in absorption coefficient (ref. 5). It can be shown that the absorption coefficient 
a is related to the light intensity by 

i = I„ Li_:«)!sa( . -. ga ( 7 ) 

1 - R 2 exp(-2o!t) 


where I is the intensity of the beam after passing through the sample, Iq is the inci- 
dent beam intensity, R is the reflectivity of the sample surface, and t is the sample 
thickness. The transmission r is the ratio of the intensities 

T _ JL _ (1 - R) 2 exp(-o?t) (g- 

IQ 1 - R 2 exp(-2a;t) 


Solving equation (8) for a gives 


a = In 


2-d2 


(1 - R) 2 + j{l - R) 4 + 4t*R 
2 r 


.-1 


( 9 ) 


It is readily seen from equation (9) that if R « 1 , a varies as the logarithm of 
the _r_e.cip meal- of _the -transm is s ion.- ~Therefore ? -as — r— app roache s-zer o ; it can be com e — 
very important that measurements of t become more accurate; yet, the sensitivity of 
the experimental setup is such that more accurate values of percent transmission cannot 
be measured. Since the present study deals primarily with energy locations and relative 
absorption strengths, values of absorption coefficients would not lead to any significant 
new conclusions and percent -transmission measurements would be just as informative. 
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Electrical Measurements 


Hall samples were cut in the standard bridge shape shown in figure 3. Contacts to 
the samples were made by sputtering aluminum on the tabs and annealing at 600 u C for 
15 minutes to form an eutectic bond between the aluminum and silicon. The leads were 


1. 3 cm 



Figure 3 • - Typical Hall shape with dimensions. 

soldered to the aluminum -coated tabs with indium solder and a tesla coil was arced to 
each tab to help break down any rectifying contacts that still remained. Resistivity and 
Hall coefficient were determined in a manner essentially the same as that described in 
reference 19. Measurements were taken reversing both current and magnetic field to 
reduce error. The magnetic field used for the Hall measurements was 2740 G. 

The resistivity of the optical samples prior to polishing was measured by a four- 
point probe as explained by Valdes in reference 38. The resistivity of some optical 
samples was not measured due to difficulty in establishing good electrical contacts. 

Doping Procedures 

Beryllium.- Beryllium was thermally diffused into the silicon using the so-called 
sandwich method of reference 30. A layer of beryllium was vapor deposited onto each 
silicon wafer. Two or more wafers were stacked with the plated sides adjacent and 
placed in a helium environment at temperatures of 1300° C to 1340° C for 1 to 2 hr (see 
table 1). The wafers were welded together during the doping procedure and had to be 
sawed apart. 

Lithium.- The most efficient way to introduce the lithium was to plate a so-called 
dummy wafer on both sides with lithium and sandwich this dummy between two beryllium- 
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doped wafers. The sandwich was placed in a helium environment and heated at 600° C 
for 30 min. The two outer wafers were then turned over and heated for another 30 min. 
The lithium -plated wafer was then removed, and the two doped wafers were annealed at 
600° C for 30 min to assure a more uniform doping. 

Annealing and Quenching Procedures 

Samples to be annealed were placed between two high-purity wafers of silicon in a 
helium environment. This was done to reduce surface chemical reactions between stray 
impurities in the system and the doped sample. 

For quenching, the samples were held against a cylindrical piece of graphite by 
means of a reduced pressure inside the graphite cylinder (see fig. 4). A thermocouple 
was placed against the silicon to allow temperature monitoring. This system was placed 
in a cylindrical quartz tube with a nitrogen gas stream flowing by the sample. The graph- 
ite was heated by means of an rf heater until the wafer reached the desired temperature 
where it was held for about 3 to 5 min and dropped into a liquid nitrogen bath. In all cases 
the samples were brought back to room temperature before mounting in the cryostat. 
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RESULTS AND DISCUSSION 


Beryllium 

When beryllium is thermally diffused into silicon, the resultant p-type material is a 
consequence of acceptor levels having an ionization energy of about 170 meV as determined 
from Hall effect measurements. The parameter plotted in figure 5 is derived from Hall 



Figure 5.- Parameter containing carrier concentration and temperature 
as a function of temperature for beryllium-doped silicon. 


coefficient calculations, and the slope of the curve in the intermediate temperature region 
is directly proportional to the ionization energy. Hall mobility, 
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(10) 


.. _ r h 
mh- — 

as a function of temperature typically behaves as shown in figure 6. 



Temperature, K 

Figure 6. - Hall mobility as a function of temperature 
for beryllium-doped silicon. 

The temperature dependence of the mobility indicates the dominant free -carrier 
scattering mechanism. If the scattering is primarily due to phonons, the mobility has a 
negative power dependency on the temperature; however, if ionized impurity scattering 
is dominant, the mobility will be a positive power function of the temperature (ref. 39). 
Here it is noted that the mobility is primarily controlled by phonon scattering down to 
temperatures around 80 K. When silicon is heavily compensated (on the order of 
10 16 cm -3), the mobility is primarily controlled by phonon scattering down to tempera- 
tures _of_only_ 100_K_or^greater;_the ref ore, -the -beryllium-is-typical-of a-p -type -impurity -- 
with low or moderately low compensation levels. 


Absorption studies in the infrared region of beryllium in silicon at low tempera- 
tures show that there are two series of sharp absorptions, located around 1450 cm - l and 
1080 cm _ l (labeled Set I and Set II, respectively) and a single broad absorption centered 
on 500 cm - l (see fig. 7). The spacings between absorptions are typical of the spacings 
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Wavelength, pm 



1500 1400 1300 1200 1100 1000 v 600 500 400 

Wave number, cm " 1 


Figure 7.- Typical beryllium spectra showing Set I and Set II absorptions 
and the broad absorption band at 500 cm" 1 . 

observed for shallow acceptors in silicon as seen in table 2 and figure 2, If shallow- 
acceptor theory is applied, these two series correspond to ionization energies of 191 
and 145 meV. This is in acceptable agreement with the value obtained from the Hall 
measurements. 

As the temperature of the sample is raised, the Set I and Set n lines broaden quite 
rapidly as shown in figures 8(a) and 8(b). Resolution of the lines above 100 K was not 

TABLE 2.- ENERGY LOCATION AND SPACINGS OF LINES IN TYPICAL 
BERYLLIUM AND BERYLLIUM-LITHIUM SPECTRA 
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Wavelength, |jm 



(a) Beryllium, Set I. 


Wavelength, gm 



(b) Beryllium, Set II. 

Figure 8.- Temperature dependence of absorptions. 

possible. However, the broad line at 500 cm - * does not seem to change very much with 
temperature, indicating that it is probably a localized impurity vibrational absorption 
associated with one of the beryllium centers (refs. 24 to 26). 

As-seen-in-tabie-2-r t-he-spacings-of-the-major-lines* in the'Set randSetTI~absorptions~ 

are quite similar. The Set II series is always weaker than the Set I, never having an 
absorption coefficient greater than about one -half that of Set I. When a beryllium -doped 
sample is quenched from a temperature greater than about 600° C, the Set II series is 
drastically reduced or disappears while the Set I and 500 cm"* lines increase slightly 
(fig. 9). When a sample is slowly cooled (on the order of 2 hr) from 600° C to room 
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Wavelength, uni 



Figure 9-- Effects of quenching and annealing on beryllium. Set I and Set II. 

temperature, the Set II series absorption is increased to its maximum value relative to 
Set I. Set I and the 500 cm -1 lines are slightly decreased in strength (fig. 9). 

From these data it appears that the Set- 1 and Set II series arise from two differ- 
ent configurations of the beryllium impurity atoms. Since the spacings between the lines 
are essentially the same and the ratio of the absorption coefficients can be changed, it 
is very unlikely that the Set II absorptions are the P-y /2 s P ec t ra as might be suspected 
at first. That is, transitions have previously been observed to states associated with 
the P 3/2 valence band and the P\f 2 valence band which are separated by about 44 meV 
(~355 cm-1) (see fig. 1). However, it is unlikely that the spacings between lines would be 
so close or that the relative strengths of absorption would be affected by heat treatment. 

It might be noted that beryllium is not unique in introducing more than one indepen- 
dent impurity level spectrum. Sulfur, when introduced into silicon, gives rise to four 
donor levels, at least two of which are unrelated (i.e., not the second level of a helium 
model), and whose relative absorptions are dependent on temperature histories (ref. 35). 

Figure 9 indicates that the 500 cm - * band is associated with the Set I impurity 
vibration as opposed to Set II since it generally tends to follow the behavior shown by 
Set I rather than Set n under quenching and annealing. When a beryllium -doped sample 
is isothermally annealed at temperatures greater than 1000° C, the Set I and 500 cm -1 
absorptions decrease at about the same rate, so that after about 4 to 6 hr they both com- 
pletely disappear due to out diffusion of the beryllium. 
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When beryllium is introduced into Czochralski-grown silicon (oxygen content typi- 
cally 1()18 cm-3) 5 the 9-jum absorption associated with interstitial oxygen completely 
disappears (fig. 10). The beryllium spectrum appears to be no different from that for the 
float-zone silicon (oxygen content less than 10 16 cm -3). Even after annealing for several 
hours to remove all the optically active beryllium, the 9-jum oxygen line is never fully 
recovered (i.e., the interstitial oxygen content in a position giving rise to the 9-pm line 
remains less than 0.1 its original value). 


Wavelength, pm 



Figure 10.- Effect of beryllium on oxygen 
in Czochralski-grown silicon. 


Beryllium -Lithium Complexes 

When lithium is thermally diffused into silicon, it normally occupies an interstitial 
position and gives rise to a donor level. It is expected that lithium would compensate the 
beryllium impurities and allow study of the singly ionized beryllium states, expected if 
beryllium could be described by a helium model. After the lithium was introduced into 
beryllium -doped silicon, the samples were still p-type after diffusion for 1 hr at 600° C. 
Calculations from Hall data (fig. 11) indicated that the holes were bound in a level with an 

ionization -energy— E-A-—oLabout—108-meV^ — This- may be- compared with- ionization-energies-- 

which would be expected from compensation of the helium model, namely energies around 
800 meV (=*6450 cm _ l). The mobility data shown in figure 12 exhibited very little of the 
ionized impurity scattering which would be expected for a highly compensated sample; in 
fact, the mobility in beryllium -lithium samples was essentially the same as for compar- 
ably doped beryllium samples. 
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Figure 11.- Parameter containing carrier concentration and temperature 
as a function of temperature for a beryllium-lithium sample. 

The suspicions that the introduction of lithium had caused something unusual were 
verified when the infrared spectrum of a beryllium-lithium-doped sample was observed. 

As shown in figure 13, two new series of acceptorlike absorptions were observed at ioni- 
zation energies corresponding to 106 and 81 meV, respectively. The previously observed 
beryllium spectra were still there. In addition, the absorption around 500 cm - * was 
asymmetric and had its minimum value at about 516 cm"^. 

The spacings between the lines of all four excited-state series are essentially the 
same (cf. table 2). After quenching a beryllium -lithium series from 600° C or greater, 
the intensities of the Set II lines are greatly reduced while the intensities of the beryllium - 
lithium series, Set I increase a little (see fig. 14). Similarly with a slow cool from 
600° C, the beryllium -lithium absorptions, Set II are increased while the beryllium - 
lithium absorptions, Set I decrease very little. 
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Percent transmission 



Temperature, K ... . i 

Wave number, cm 1 


Figure 12.- Hall mobility as a Figure 13 . - Typical spectra showing beryllium- 

function of temperature for lithium absorptions, Set I and Set II. 

a beryllium-lithium sample. 


600°C quench 

600°C slow cool 



Wave number, cm‘1 


Figure 14. - Effects of quenching and annealing on beryllium-lithium absorptions. 

Set I and Set II. 
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As with the beryllium samples, the absorptions associated with beryllium -lithium 
broaden rapidly with increasing temperature and are unresolvable at temperatures 
greater than about 60 K to 70 K. However, as the sample temperature is increased, 
another series of lines appear which are displaced from the beryllium -lithium series, 
Set I by about 1.2 meV (~10 cm _ l) (see fig. 15). 


Percent transmission 
(arbitrary scale) 



Wave number, cm' 1 

Figure 15*- Temperature dependence of beryllium-lithium absorptions. 

Set I and Set II. 

The absorptions associated with beryllium -lithium, Set II do not seem to broaden 
as rapidly and do not show the second series. The appearance of the small absorptions 
in the beryllium -lithium, Set I is consistent with what would be expected if the fourfold 
degenerate ground state were split into two levels separated by 1.2 meV. As the tem- 
perature is raised, the upper level is populated by thermal excitation of the holes, and 
the transitions to the excited states can then be observed. 

When lithium is introduced into a beryllium -doped Czochralski -silicon sample, the 
beryllium -lithium spectra are essentially the same as for the float -zone material. The 
spectra for a Czochralski sample are shown in figure 13. The beryllium -lithium, Set II 
lines are broad and lie in the vicinity of strong lattice absorptions which make accurate 
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determination of the spacings very difficult. It is felt that comparative spectra may be 
the best way to study accurately this series of absorptions. 

PROPOSED MODELS 
Beryllium 

The very fast diffusion rate of beryllium in silicon indicates interstitial diffusion 
and leads one to expect beryllium to be an interstitial impurity. However, the addition 
of beryllium produces p-type silicon which indicates that the electrically active beryl- 
lium is in a substitutional position. Mass -spectroscopy data show that only about 10 per- 
cent of the beryllium contributes to the change in electrical and optical properties. It is 
believed that the other 90 percent of the beryllium is located in neutral sites, most likely 
close substitutional -interstitial beryllium pairs which would behave as neutral impurities. 
It appears that the beryllium diffuses interstitially but is trapped in the vacancies which 
are abundant at the doping temperature of 1300° C. When the beryllium is introduced 
into Czochralski-grown silicon, the interstitial beryllium interacts with the interstitial 
oxygen to form an interstitial beryllium -oxygen molecule. The bond energy of the 
beryllium -oxygen and the silicon-silicon bonds exceeds that of the silicon-oxygen-silicon 
bqnds making the formation of the interstitial beryllium -oxygen energetically favorable. 

Set I.- If a beryllium impurity atom were located on a substitutional lattice site, 
it would satisfy two of the covalent bonds with nearest neighbor silicon atoms. The 
remaining covalent bonds would capture electrons, but the entire center would remain 
neutral by binding two holes to the beryllium center on the lattice site which has a double 
negative charge (fig. 16). (Figures 16 to 21 are oversimplified for pedagogical reasons. 

In fact, the trapped electrons cannot be localized on any particular bond and the holes 
would be in orbit around the impurity center.) The holes would be several atomic diam- 
eters away and would be described in a simple way by effective -mass theory for a helium 
model. Therefore, when one of the holes is raised to an excited state, the spectra should 
look essentially like a Group HI shallow acceptor. This is presumably what gives rise 
to the Set I absorptions. In the present context, the beryllium atom appears to be sig- 
nificantly smaller than the silicon atom. According to the theory of Morgan (ref. 15), 
this smaller atom will have a large local strain field associated with it which produces 
ground stat es muc h deeper Jhan predicted by -effective -mass theory . ' 

When one of the holes is freed, leaving an ionized center, simple theory would pre- 
dict that the second hole should be observable at an energy somewhere around four times 
that required to free the first hole. This type of behavior has been observed for the 
donors magnesium and sulfur in silicon and the acceptor zinc in germanium. In the case 
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Figure l6.- Schematic model for neutral beryllium. Set I. 


of beryllium, however, no such second level has been observed in this region. Several 
things could explain the absence of the second series of absorptions associated with the 
Set I lines. 

One might start over by assuming that the beryllium is not in a substitutional site. 
This would require entirely new model concepts for the beryllium acceptor and would 
presumably complicate the model for the beryllium -lithium complex. On this basis, 
another explanation was sought. When the first hole is freed from the heliumlike model, 
the second hole will be attracted into a smaller orbit (see fig. 17). It is clearly possible 
that when the orbit becomes smaller, the assumption of a dielectric constant used in 
effective -mass theory becomes much weaker. The screening now seen by the hole could 
be much less than before, resulting in a much more tightly bound hole than predicted. 

In addition, a lattice -symmetrical short-range potential (central-cell potential) 
would be much more effective for a hole located close to the impurity site. This potential 
would be used for calculating the ground-state energy for various impurities, and large 
values of the potential might possibly result in more tightly bound holes. In a theoretical 
paper, Ryabokon and Svidzinskii (ref. 14) state that for Group II impurities, ". . . the 
effective short-range potential turns out to be more effective in them and leads to quite 
large level shifts." Furthermore, they add that some of the deep impurity levels in sili- 
con fall inside allowed energy bands. Thus, it seems possible for beryllium to be in a 
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Figure 17. - Schematic model for singly ionized beryllium. Set I. 

substitutional site, even though the second hole is not observed as predicted by the helium 
model. 

Set II .- When a beryllium sample is quenched from high temperatures, the Set II 
absorptions are greatly reduced or disappear, and the Set I absorptions are enhanced. 

After a slow cool, the Set II absorptions are much stronger, and the Set I absorptions 
are slightly reduced. This seems to indicate that the Set II absorptions result from a 
more complex arrangement than Set I and that Set I absorptions are due to a component 
of impurities giving rise to Set II. It appears that the simplest model for the Set II 
absorptions is two beryllium atoms on nearest neighbor lattice sites. This configuration 
is somewhat similar to the divacancy bonding configurations, but less complex (see fig. 18). 
As the temperature of the sample is raised, the thermal energy of the lattice is such that 
the Set II sites are broken up into two Set I sites, that is, two substitutional beryllium 
sites. This results in a reduction of Set II absorptions after a quench. When the sample 
is slo wly co.ole.d.-the-bery.llium-atoms-pr-esumablv-tend-to-aggregate~to~relieve~crvstal 
stresses around the small beryllium impurity. Thus, more pairs are formed, and Set II 
absorptions are stronger after a slow cool. Again, the second hole associated with this 
model might not be observable in the infrared range (see fig. 19). These two models 
seem to be consistent with the optical and electrical data. 
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Beryllium and Lithium 

Set I.- When lithium is introduced into beryllium -doped silicon, an acceptor level 
with essentially the same spacings as the beryllium, Set I spectrum has been observed, 
but with the binding energy of the hole reduced by about 45 percent. It is believed that 
the lithium goes into one of the nearest neighbor interstitial tetrahedral sites (ref. 21) 
and its valence electron satisfies one of the unsatisfied covalent silicon bonds (fig. 20). 

The lithium in this position could reduce the local strain field, thus resulting in a lower 
ionization energy (ref. 15). While the short-range effects of the lithium are apparently 
large, the long-range effects should be quite small; therefore, the excited-state spectra 
should be essentially unaffected. Similarly, if the lithium is located in the immediate 
vicinity of the beryllium impurity, the ionic scattering from the complex center would be 
less than if the spacing between the two centers were large. That is, the beryllium - 
lithium center would be essentially a monovalent acceptor site. 

As pointed out earlier, there is an apparent splitting of the ground state for the 
loosely bound hole in beryllium -lithium samples. If the lithium is in the immediate vicin- 
ity of the beryllium, the interaction between the beryllium and lithium ions would most 
likely be directed along a bond or <111> axis. In effect, this would be equivalent to a uni- 
axial stress applied in the <111> direction. As reported in reference 8, a <111> stress 
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Figure 20.- Schematic model for neutral beryllium-lithium complex, Set I. 


26 





























Figure 21.- Schematic model for neutral beryllium-lithium complex, Set II. 


causes the ground state to split into two levels. Since these stresses are local, the 
excited states would not be split. This is in agreement with what is observed. 

Set II.- When lithium comes into the vicinity of the beryllium, Set II, its valence 
electron again goes into one of the unsatisfied silicon convalent bonds, removing what 
would correspond to the deep hole (see fig. 21). The result, as before, would be to reduce 
the binding energy of the remaining hole. The mobility would not be affected by this con- 
figuration because no new ionic scattering centers have been introduced. 

If this model is accurate, the complex would probably not have a <11 1> symmetry 
when the lithium is introduced; but the model for the beryllium complex, Set II would have 
a <11 1> axis of symmetry. However, the ground state for a heliumlike model would not 
have the same degeneracy as the Group HI acceptors; and in fact, it is likely that both the 
beryllium, Set I and Set II will have only twofold degenerate ground states. The lithium 
should reduce even further the symmetry of the beryllium, Set II center, but the ground 
state would probably remain twofold degenerate (ref. 8). 


CONCLUDING REMARKS 


Beryllium as an impurity in silicon gives rise to two acceptor levels and associated 
excited-state levels referred to as Set I and Set II. The ratio of the absorption coeffi- 
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cients between Set I and Set II can be changed by quenching or annealing. Analysis of 
these studies indicates that the Set II absorptions are due to an impurity complex and that 
the Set I absorptions are due to a component of the Set II complex. 

When lithium is introduced into the beryllium -silicon system, it gives rise to two 
new acceptor levels. These levels behave under quenching and annealing just like the 
beryllium, Set I and Set II levels. The beryllium -lithium, Set I ground state is split into 
two levels. 

The following models, which seem to be the simplest impurity configurations that 
agree qualitatively with all the observations, are proposed to explain these four defects: 

1. For beryllium, Set I, this defect is a single substitutional beryllium impurity. 
The second hole of a heliumlike model is so tightly bound that it cannot be found in the 
region predicted by shallow -acceptor theory. 

:) 2. For beryllium, Set II, this defect is made up of two beryllium atoms on nearest 

neighbor substitutional sites. Again, a deep hole may be associated with the level; but 
the hole could be so tightly bound that it falls inside the allowed energy bands. 

3. For beryllium -lithium, Set I, the lithium ion is trapped in the immediate vicinity 
of the beryllium, Set I and removes the deep hole. The symmetry of the impurity site is 
reduced, resulting in the reduction of the degeneracy of the ground state and the splitting 
of this level. 

4. For beryllium -lithium, Set II, the lithium removes the deep hole as in Set I; but 
the ground-state degeneracy is not reduced because of the different symmetries asso- 
ciated with the beryllium, Set I and Set II sites. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., November 20, 1972. 
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